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Abstract To understand the neural origins of rhythmic behavior one must characterize the 
central pattern generator circuit and quantify the population size needed to sustain functionality. 
Breathing-related interneurons of the brainstem pre-Botzinger complex (preBotC) that putatively 
comprise the core respiratory rhythm generator in mammals are derived from Dbx 7-expressing 
precursors. Here, we show that selective photonic destruction of Dbx1 preBotC neurons in neonatal 
mouse slices impairs respiratory rhythm but surprisingly also the magnitude of motor output; 
respiratory hypoglossal nerve discharge decreased and its frequency steadily diminished until 
rhythm stopped irreversibly after 85±20 (mean ± SEM) cellular ablations, which corresponds to -15% 
of the estimated population. These results demonstrate that a single canonical interneuron class 
generates respiratory rhythm and contributes in a premotor capacity, whereas these functions are 
normally attributed to discrete populations. We also establish quantitative cellular parameters that 
govern network viability, which may have ramifications for respiratory pathology in disease states. 
DOI: 10.7554/eLife.03427.001 



Introduction 

Central pattern generator (CPG) circuits give rise to common behaviors such as swimming, walking, 
and breathing {Grillner, 2006; Grillner and Jessell, 2009; Kiehn, 2011). To understand the cellular 
origins of these behaviors, key problems are to identify the rhythmogenic and premotor populations, 
and then quantify the requisite number of neurons to sustain functionality. We address these issues 
in the mammalian breathing CPG by cumulatively ablating a genetically identified interneuron pop- 
ulation hypothesized to form the rhythmogenic core while monitoring effects on network output in 
real time. 

The brainstem pre-Botzinger complex (preBotC) putatively drives inspiratory breathing rhythms 
(Smith etaL, 1991; Feldman etaL, 2013; Moore et al., 2013). These rhythms persist in reduced slice 
preparations that retain the preBotC and can be monitored via respiratory hypoglossal (XII) nerve dis- 
charge, providing a powerful in vitro model of breathing behavior {Lieske et al., 2000; Koizumi et al., 
2008; Funic and Greer, 2013). Rhythmogenic preBotC interneurons are distinguished by glutamater- 
gic transmitter phenotype or the expression of neuropeptides and peptide receptors, or their inter- 
section [Funk et al., 1993; Gray et al., 1999, 2001; Guyenet et al., 2002; Stornetta et al., 2003; 
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eLlfe digest Our first breath, moments after we are born, is the result of a pattern of activity in 
our brain that started in the embryo and will continue almost effortlessly until we die. Like other 
rhythmic activities, such as walking and swimming, breathing originates from circuits of neurons in 
the brain that generate patterns. These circuits pass messages to other cells that translate them 
into the physical movements required to take a breath. Interrupting these patterns by injury or 
illness can lead to breathing disorders or cause death. 

Previous studies have identified a class of neuron, which all express a specific gene, that is 
necessary for breathing. Mice born without this class of cell failed to ever take a breath and died at 
birth. These neurons are found in part of the brainstem and can continue to generate rhythm even 
when this section of the brainstem is removed from newborn mice and cut into very thin slices. 
However, it is unclear how many of these neurons are needed to maintain a breathing rhythm. 

Wang et al. used a laser to destroy the breathing rhythm-generating neurons in these slices one 
at a time and found that the rhythm of breathing in (i.e., inspiration) stopped after -15% of the neurons 
were destroyed. This suggests that a high percentage of these neurons must be maintained for 
breathing to continue normally. 

Wang et al. also discovered that destroying the rhythm-generating neurons reduced the strength 
of the signals sent from the brainstem to trigger the movements that cause breathing in. This 
suggests that the same class of neurons also sends messages to the muscles involved in breathing; 
it was previously thought that a separate class of cell in the same part of the brain sent these 
messages. 

Studies involving live animals are now needed to confirm the results. If confirmed, the findings 
may be used to develop new treatments for a number of breathing disorders. Medications that 
boost the signals sent to the muscles by these neurons might be useful for treating sleep apnea. 
Wang et al. also suggest that medications that boost rhythm generation might be useful for 
premature infants with breathing difficulties and people with drug-induced breathing problems. 
Moreover, finding ways to maintain breathing rhythms with fewer of these neurons may help those 
with neurodegenerative disorders, which cause cells in the brain to be lost. 
DOI: 10.7554/eLife.03427.002 



Wallen-Mackenzie et al., 2006; Tan eta/., 2008). Glutamatergic, peptidergic, and peptide receptor- 
expressing preBotC interneurons develop from a common set of precursors that express the home- 
obox gene Dbx1 during embryonic development. Dbx7-derived interneurons (hereafter Dbx1 neurons) 
in perinatal mice are inspiratory modulated, and Dbx7-null mice die at birth without ever breathing 
{Bouvier et al., 2010; Gray et al., 2010; Picardo et al., 2013). Therefore, we — and others — proposed 
that Dbx1 preBotC neurons comprise the core inspiratory rhythm generator, i.e., the Dbx1 core 
hypothesis. 

Previously, we estimated the cellular mass critical for respiratory rhythm generation by laser- 
ablating preBotC inspiratory interneurons identified by Ca^^ imaging. The destruction of -120 
rhythmic neurons irreversibly stopped respiratory rhythmogenesis {Hayes et al., 2012). However, 
inspiratory-modulated preBotC neurons may be excitatory or inhibitory. Ca^^ fluorescence changes 
cannot differentiate rhythmogenic glutamatergic neurons {Funk et al., 1993; Wallen-Mackenzie 
et al., 2006) from GABA- or glycinergic neurons {Kuwana et al., 2006; Winter et al., 2009), which 
influence sensory integration and coordinated patterns of muscle contraction for respiratory behav- 
iors, but are non-rhythmogenic because inspiratory rhythms in vivo and in vitro do not require synaptic 
inhibition {Feldman and Smith, 1989; Brockhaus and Ballanyi, 1998; Kuwana et al., 2003; Ren and 
Greer, 2006; Fujii et al., 2007; Janczewski et al., 2013). Therefore, we predicted that the selective 
destruction of Dbxl preBotC neurons, which are predominantly glutamatergic {Bouvier et al., 2010; 
Gray et al., 2010) — unlike locomotor Dbxl neurons in lumbar spinal cord of which -70% express 
inhibitory transmitters {Lanuza et al., 2004; Talpalar et al., 2013) — would impair rhythmogenesis with 
a cell-ablation tally much lower than 120 {Hayes et al., 2012). To test this prediction of the Dbxl core 
hypothesis, we used photonics to map the positions of Dbxl preBotC neurons, and then laser ablated 
them — one cell at a time — while continuously monitoring respiratory motor output {Wang et al., 2013). 
As predicted, cumulative destruction of Dbxl preBotC neurons progressively decreased respiratory 
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frequency until rhythm ceased after ablation of -85 neurons. Surprisingly, cumulative Dbx1 cellular 
ablations also diminished the amplitude of respiratory XII nerve discharge, suggesting that Dbx1 
preBotC neurons also influence motor output. In simulations that assign only rhythm-generating func- 
tion to preBotC neurons, cumulative ablations decreased frequency and stopped rhythmogenesis but 
at much lower tallies and without perturbing the output amplitude. These ablation results and model- 
experiment discrepancies, combined with antidromic activation from the XII nucleus and axonal pro- 
jection patterns, ascribe rhythm-generating and premotor roles to Dbx1 preBotC neurons. Thus, we 
demonstrate that one cardinal class of hindbrain interneurons serves two distinct roles in a key mam- 
malian CPG; the data further establish quantitative cellular parameters that minimally ensure network 
functionality. 

Results 

Ablation of Dbx1 preBotC neurons precludes rhythmogenesis 

Dbx1 neurons were detected and mapped within the preBotC, and then laser ablated individually, in 
sequence, while monitoring respiratory network functionality via XII motor nerve output. Experiments 
began with an 'initialization phase' that defined the domain for detection and ablation, which 
was bilateral. We used a Dbx1 Cre-driver line (Dbx7*^'^^'^^^) coupled with floxed reporter mice 
{Gt(ROSA)26Sor^'°'''^^°P'^^^°'^^^°) to locate Dbx1 neurons via fluorescence. Viewed in the transverse plane 
of slices that expose the preBotC at their surface (i.e., preBotC-surface slices), Dbx1 neurons form a 
bilaterally symmetrical V-shape starting dorsally at the border of the XII motor nuclei and continuing 
ventrolaterally to the preBotC {Figure 1A and Figure 1 — figure supplement 1A). Dorsally, the preBotC 
adjoins the semi-compact division of the nucleus ambiguus (scNA); the ventral border of the preBotC 
is orthogonal to the dorsal boundary of the principal sub-nucleus of the inferior olive (lOPioop) 
{Ruanglcittisalcul et al., 2011, 2014). These spatial relationships visible in bright field or epifluorescence 
allow us to pinpoint the preBotC {Figure 1A). At the cellular level, identifying putative rhythmogenic 
neurons on the basis of fluorescent protein expression alone is acceptable because the overwhelming 
majority of Dbx1 preBotC neurons are inspiratory (e.g.. Figure IB) {Picardo et a!., 2013). 

In the subsequent 'detection phase', a visible wavelength laser scanned the domain and an iterative 
threshold-crossing algorithm analyzed the image to then draw regions of interest (ROIs) for putative 
cell targets based on fluorescence brightness. Potential targets were evaluated on the basis of shape 
to differentiate cell bodies from auto-fluorescent debris, and to reject the fluorescence from dendrites 
and neuropil whose somata were detectable in adjacent focal planes (Wang et al., 2013) {Figure 1 — 
figure supplement 2). The map of ROIs for validated cell targets was retained at each focal plane 
{Figure 1C — figure supplements IB and 2, red ROIs). Potential targets that did not meet these 
criteria were discarded (but displayed for demonstration purposes in Figure 1C and Figure 1 — figure 
supplements IB and 2, blue ROIs). Target detection was repeated at lO-pm increments through the 
z-axis and the final three-dimensional map of targets was stored in memory {Figure ID). Typically, we 
detected 26-50 Dbx1 neurons per focal plane per side {Figure 1 — figure supplement 3) for a total 
average number of 705 targets in the preBotC (SD 1 1 9, SEM 59, range: 548 to 802, n = 8 slices). 

During the 'ablation phase' of the experiments, Dbx1 neurons in preBotC-surface slices were ran- 
domly selected for photonic lesioning. Each target in the domain was individually spot scanned with a 
Tiisapphire laser using maximum intensity 800-nm pulses until target destruction was confirmed by 
three forms of optical criteria (Wang et al., 2013) or was deemed a failure. Generally >90% of lesion 
attempts are successful {Figure 2 — figure supplement 1) {Hayes et al., 2012; Wang et al., 2013). 
Only confirmed lesions add to the running tally. 

The frequency and amplitude of inspiratory motor output diminished at the onset of the ablation 
phase {Figures 2A and 3A). XII amplitude decreased steeply with the tally of ablated cells, and then 
stabilized at 44% of its pre-lesion value (SD 4%, SEM 1%, suction electrode recordings are reported in 
normalized arbitrary units). Frequency, however, continued to decrease (i.e., cycle period increased) 
throughout the ablation phase. Initially, within the first dozen ablations, the average decrease in respi- 
ratory frequency was nearly twofold, and it continued to fall steeply until rhythm cessation (range 
of frequencies: 0.22-0.007 Hz, Figure 3B [inset shows bi-exponential increase in cycle period], n = 5 
slices). Furthermore, the rhythm destabilized during the ablation phase. We defined regularity score 
(RS) as the ratio of the present cycle period with respect to the mean period over 10 prior cycles (see 
'Materials and methods' for RS formula). Cycle-to-cycle variations in RS indicate irregularity; the system 
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Figure 1 . Dbx1 preBotC neurons. (A) Bright field (left) and fluorescence (right) images of the right half of a preBotC- 
surface slice preparation. Anatomical landmarks are illustrated including: XII, the hypoglossal motor nucleus; scNA, 
semi-compact nucleus ambiguus; lOPioop, the dorsal loop of the principal inferior olive; and the ventral border of 
the preBotC, which is orthogonal to the lOPioop- Scale bar is 300 pm. At right, the larger white box shows the 

Figure 1. Continued on next page 
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Figure 1. Continued 

detection and ablation domain. (B) Expansion of smaller white box in A, showing tdTomato expression in Dbxl 
neurons and intracellular dialysis via patch pipette with Alexa 488 from the recorded neuron whose robust 
inspiratory discharge is illustrated at right (scale bar is 10 pm). Respiratory motor output from the XII nerve is shown 
in raw and RMS-smoothed form. Voltage and time calibration bars represent 20 mV and 2 s. Baseline membrane 
potential in the recorded neuron was -60 mV. (C) Mask of targets showing validated Dbx1 interneuron targets (red) 
and regions of fluorescence that do not pass muster and were rejected as targets (blue) for focal planes at depths 
z = (30-60 [jm). The region shown in each case maps to the 41 2 x 41 2 pm^ square shown by the larger white box 
in A (right). Only a subset of the masks are shown for economy of display. (D) 3D reconstruction of detected targets 
for all focal planes z = (10-80 pm) from the left and right preBotC. Each Dbxl neuron is represented by a single red 
point centered on its soma. 
DOI: 10.7554/eLife.03427.003 

The following figure supplements are available for figure 1 : 
Figure supplement 1. Detection of Dbxl preBotC neurons. 
DOI: 10.7554/eLife.03427.004 

Figure supplement 2. Detection of Dbxl neuron targets via fluorescence and image processing. 
DOI: 10.7554/eLife.03427.005 

Figure supplement 3. Average number of Dbxl neurons detected at each acquisition depth from z = 0 (surface) to 
z = -80 [jm in preBotC-surface slices and control slices with the ventral respiratory column (VRC) exposed at the 
slice surface. 

DOI: 10.7554/eLife.03427.006 



is trending slower when RS exceeds unity. The RS of preBotC-surface slices measured 1-9 during the 
ablation phase {Figure 3Q. Respiratory rhythm ceased altogether after an average of 85 confirmed 
Dbx1 neuron ablations in preBotC-surface slices (SD 44, SEM 20, range 42-137, n = 5 slices), well 
before exhausting the average list of 705 targets per slice. These ablations were bilateral and the tally 
reflects the sum of both sides (Figure 2 — figure supplement 2). The representative experiment in 
Figure 2A shows rhythm cessation after 62 confirmed ablations, corresponding to 9% of the total 677 
detected targets. 

Ablation of Dbx1 neurons from the ventral respiratory column 

Detection and ablation were similarly performed bilaterally in control slices whose rostral surface 
exposed the ventral respiratory column, which occupies a comparable domain for detection and abla- 
tion in the transverse plane, but this domain is -100 |jm rostral to preBotC. The ventral respiratory 
column contains inspiratory and expiratory-modulated neurons that are not associated with rhythmo- 
genesis {Smith et al., 1990; Feldman et al., 2013). 

During the initialization phase in control slices, the domain was centered on the highest density of 
fluorescent Dbxl neurons bounded by the compact division of the nucleus ambiguus (cNA) dorsally 
and the ventral margin of the slice {Figure 4). We acquired 38-60 Dbxl neuron targets per focal plane 
per side {Figure 1 — figure supplement 3) for a total average of 906 targets per slice (SD 97, SEM 34, 
range: 722-1004, n = 8 slices). 

Dbxl ventral respiratory column neurons were lesioned in random sequence during the ablation 
phase of control experiments. The amplitude of XII motor output decreased to 77% of control (SD 2%, 
SEM 1%, normalized arbitrary units) over the course of the ablation phase {Figures 2B and 3A). 
Frequency did not change. It measured 0.36 Hz (SD 0.2 Hz, SEM 0.04 Hz) during the detection phase 
compared to 0.39 Hz (SD 0.3 Hz, SEM 0.01 Hz) during the ablation phase, which was not significant 
(p=0.49, Mann-Whitney (J-test, Figure 3B). RS did not deviate from ~1 throughout the ablation phase 
{Figure 3D, n = 8 slices). These data indicate that cumulative sequential ablation of Dbxl ventral res- 
piratory column neurons has no effect on the stability or the period of the respiratory cycle. The abla- 
tion protocol exhausted the entire set of targets in every control experiment without stopping the 
rhythm (e.g., 923 confirmed ablations in Figure 2B). Ablations in control slices were also performed 
bilaterally, where the tally reflects the sum of both sides {Figure 2 — figure supplement 2). 

Transient recovery of irregular and unsustainable rhythm 

Cumulative deletion of Dbxl preBotC neurons appeared to degrade respiratory oscillator function. 
Nonetheless, an alternative explanation could involve the loss of excitatory drive (rather than destruc- 
tion of CPG core circuitry). Dbxl preBotC neurons express neurokinin-1 peptide receptors (NKIRs) 
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Figure 2. Cumulative serial ablation of Dbx1 neurons in preBotC-surface slices (A) and control slices whose surface 
exposes the ventral respiratory column, not preBotC (B). (A and B) The x-axis is a timeline. The y-axis plots XII 
amplitude (normalized units, top) and respiratory period (bottom). The respiratory period axis is continuous 
(0-200 s) but plotted with two scales. Major ticks are separated by 10 s from 0 to 20 s (with unlabeled minor ticks at 
5 s increments), and thereafter major ticks are plotted in 1 00 s divisions from 21 to 200 s (with unlabeled minor ticks 
at 50 s increments). The discontinuity in the y-axis stops at 20 s (lower portion) and starts at 21 s (upper portion). 
There is one data point for every individual respiratory period measured. The recording in A is no longer displayed 
after 6 min of XII quiescence. Substance P (SP) injection in A is displayed at higher sweep speed in Figure 5C. The 
recording in B is no longer displayed after 90 min of continuous stable XII output following the end of the ablation 
phase. Time calibrations in A and B are shown separately. 
DOI: 10.7554/eLife.03427.007 

The following figure supplements are available for figure 2: 
Figure supplement 1. Cellular laser ablation and confirmation. 
DOI: 10.7554/eLife.03427.008 

Figure supplement 2. Cumulative tally of laser ablations for preBotC-surface slices (magenta) and control slices 
whose surface exposes the ventral respiratory column, not preBotC (cyan). 
DOI: 10.7554/eLife.03427.009 



(Bouvier et aL, 2010; Gray et a\., 2010) that stimulate respiratory rhythmogenesis {Gray et a\., 1999; 
Pagliardini et aL, 2005; Ballanyi and Ruangkittisakul, 2009). Monoaminergic and peptidergic raphe 
neurons project to, and receive feedback projections from, the preBotC to elevate excitability in the 
respiratory network {Ptak et aL, 2009). Therefore, laser ablation in the preBotC could break the link 
with the raphe and thus diminish excitatory drive. To test that idea, we exposed the lesioned preBotC 
to a bolus of neuropeptide substance P (SP, 1 mM) after the respiratory cycle period exceeded 120 s. 
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Figure 3. Ablation effects on respiratory frequency and the amplitude of XII motor output. (A-D) Measurements 
are displayed in light grey and red for preBotC-surface slices and dark grey and blue for control slices that expose 
the ventral respiratory column (VRC). (A) XII amplitude and (B) respiratory frequency for preBotC-surface and 
control slices are plotted vs cumulative percent of total lesions during the ablation phase (bars show SD). Inset in 
B shows respiratory period in lieu of frequency (bars show SD) for preBotC-surface slices. (C and D) The regularity 
score {RS) is plotted vs cumulative percent of total lesions for preBotC-surface (C) and control slices (D). B, C, and D 
are plotted on semi-log axes. B and C are labeled with subordinate ticks at 2, 4, and 6. Tick labels are omitted from 
D because they match C exactly. B (inset) has linear axes. 
DOI: 10.7554/eLife.03427.010 



which we previously determined was a reliable benchmark of a slice that would cease rhythmic func- 
tion within 5-10 min without pharmacological intervention {Hayes et aL, 2012). 

First, as a control, we applied SP to unlesioned preBotC-surface slices, which transiently increased res- 
piratory frequency, i.e., lowered the cycle period from 4.7 s (SD 0.8 s, SEM 0.1 s) to 3.3 s (SD 0.5 s, SEM 0.1 
s, average period computed for 25 cycles), and then equilibrated in 21 min (SD 4 min, SEM 2 min, n = 4 
slices) {Figure 5A). The regularity score was ~1 throughout the bout, which is consistent with the stable 
rhythm expected in unlesioned slices {Figure 5B). Then, SP was injected into five preBotC-surface slices 
wherein the cumulative laser ablation of Dbxl neurons caused 120 s of quiescence. SP transiently revived 
respiratory rhythm; the average cycle period was 1.7 s (SD 0.4 s, SEM 0.2 s, computed for 10 cycles after 
SP bolus injection) but the cycle period slowed down rapidly, surpassing the control period previously 
measured during the detection phase (4.7 s) within 3 min (SD 2 min, SEM 1 min). Cycle period continuously 
lengthened and fluctuated from cycle to cycle, and then the rhythm stopped altogether {Figure 5Q. 
Judged on the basis of equilibration time, the transient effects of SP were significantly briefer in lesioned 
slices (p=0.02, Mann-Whitney (J-test). Furthermore, lesioned slices ultimately fell inexorably silent {Figures 
2A and 5Q, whereas unlesioned slices maintained rhythmicity for 4-6 hr {Figure 5A). More importantly, the 
SP-evoked activity in lesioned preBotC-surface slices was irregular: RS measured 2-1 0 {Figure 5D, compare 
to unlesioned slice in Figure 5B). These data indicate that NK1 R-expressing Dbxl neurons can evoke 
transient cycles of respiratory activity, but the loss of -85 Dbxl preBotC neurons slows the respiratory 
oscillator frequency and renders XII motor output nonfunctional. 
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Figure 4. Dbx1 neurons in the ventral respiratory column. (A) Bright field (left) and fluorescence (right) images of 
the right half of a control slice preparation. Anatomical landmarks are illustrated including: XII, the hypoglossal 
motor nucleus; cNA, the compact division of the nucleus ambiguus; lOPioop, the ventral portion (loop) of principal 
Figure 4. Continued on next page 
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Figure 4. Continued 

sub-nucleus of the inferior olive; and VRC, the ventral border of the ventral respiratory column. Scale bar is 300 pm. 
At right, the larger white box shows the detection and ablation domain. (B) Expansion of smaller white box in A, 
showing tdTomato expression in Dbxl ventral respiratory column neurons (scale bar is 10 pm), one of which was 
recorded. Intracellular dialysis via patch pipette with Alexa 488 is visible in the recorded neuron whose inspiratory 
depolarization and discharge pattern are illustrated at right. Respiratory motor output from the XII nerve is shown 
in raw and RMS-smoothed form. Voltage and time calibration bars represent 20 mV and 1 s. (C) Masks of targets 
showing validated Dbx1 interneuron targets (red) and regions of fluorescence that do not pass muster and were 
rejected as targets (blue) for focal planes at depths z = (40-70 pm). (D) 3D reconstruction of detected targets for all 
focal planes z = (0-80 pm) in the ventral respiratory column (VRC) of the left and right side. A single red point 
centered on its soma represents each Dbx1 neuron. The highest fraction of accepted Dbxl target cells is found at 
deeper focal planes (see Figure 1 — figure supplement 2 and 'priority rule' explained in 'Materials and methods'). 
DOI: l0.7554/eLife.03427.011 



Modeling Dbx1 neuron ablation in the preBotC 

We used graph theory and simulations to investigate how Dbx1 neuron ablations affect preBotC 
structure and function. The Rubin-Hayes preBotC neuron model {Rubin et aL, 2009) was assembled 
in Erdos-Renyi G(ap) graphs (Newman etaL, 2006) with population sizes n from 200-400 and connection 




Figure 5. Substance-P (SP) injections in preBotC-surface slices. (A) SP bolus injected in an un-lesioned preBotC- 
surface slice. XII output magnitude is plotted with cycle period as a time series. (B) Semi-log plot of regularity score 
[RS] for 30 min after SP injection from the slice preparation in A. R5 axis is continuous but plotted with two scales. 
(C) preBotC-surface slice shown in the acquisition phase (left) and during the ablation phase (right), which were 
separated by a time gap of 3 hr. After 1 20 s of quiescence (data point circled in red), SP injection revived the 
rhythm transiently. (D) Semi-log plot of RSfor 15 min after SP injection from the slice preparation in C. Data in C 
and D were from the same preparation as in Figure 2A. 
DOI: 10.7554/eLife.03427.012 



Wang et a/. eLife 2014;3:e03427. DOI: 10.7554/eLife.03427 



9 of 25 



Research article 



Neuroscience 



probabilities p from 0.1 to 0.2. These parameter ranges encompass n = 325, an empirical estimate of 
the number of excitatory neurons in the perinatal mouse preBotC (Hayes et aL, 2012) as well as 
p=0.13, the only experimentally determined connection probability among putatively rhythmogenic 
preBotC neurons in acute mouse slices {Rekling et aL, 2000). Networks within the above n-p param- 
eter range that generated respiratory-like cycle periods of ~4 s are shown with asterisks in Figure 6A 
and Figure 6 — figure supplement 1A. This set of model networks also generated network-wide 
bursts within 200-300 ms following brief focal glutamatergic stimulation of five or more constituent 
neurons {Figure 6B) in agreement with focal glutamate un-caging experiments in neonatal mouse 
slices, which showed that simultaneous stimulation of 4-9 preBotC neurons can trigger inspiratory 
bursts with similar latency {Kam et al., 2013b). These results substantiate that the model networks 
well represent the neonatal mouse preBotC in vitro. 
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Figure 6. Numerical simulations. (A) Networks of Dbxl preBotC neurons with population size (n) and synaptic connection probability (p). Blocks show 
the mean cycle period according to the colorimetric scale (right) for 10 (or more) realizations of the network for each (ap) pair. Asterisks denote networks 
that generated respiratory-like (~4 s) cycle periods in >80% of individual realizations of the network. (B) Focal glutamatergic stimulation of constituent 
neurons in a model network (n,p) = (330,0.125). Network-wide bursts can be evoked when five or more individual cells are stimulated. These simulations 
mimic holographic laser-mediated glutamate un-caging experiments [Kam et aL, 2013b) and are included because they bolster confidence that our 
model networks accurately capture features and behaviors of the preBotC in newborn mice. Raster plots show spike activity in six constituent neurons 
randomly selected from the network and focally stimulated (see 'Materials and methods' for numerical simulation of glutamate un-caging protocol). If 
focal stimulation evoked EPSPs (not spikes) then the raster reports 'EPSPs'; spikes are indicated by short vertical lines. From left to right, the number of 
stimulated units increments by one; five (or more) units evoked an inspiratory-like burst. A running-time histogram of network activity is shown at the 
bottom. Calibration bars represent 1 00 spikes/1 0-ms bin (vertical) and 0.5 s (horizontal). (C) Running-time histogram for one simulation of sequential 
ablation in a network (n,p) = (330, 0.125). Cell ablation tally is shown (top). Time calibration is 30 s. Spikes-per-bin calibration bar is the same as the inset 
(lower), 100 spikes/1 0-ms bin. Insets show a raster plot of spike activity in the entire network with a running-time histogram. The numerical y-axis reports 
cell index for each neuron model in the network. Left inset shows the first ablation (magenta arrow). Right inset shows all cumulative 36 ablations 
(magenta arrows). Time calibration for both insets is 1 s (at right). (D) Cycle period and (E) spikes-per bin (i.e., a measure of the magnitude of simulated 
network output as in C) are plotted vs cumulative percent of total ablations for 10 networks with (n,p) = (330, 0.125). D plotted in semi-log axes, E in 
linear axes. Magenta shows data from individual networks, cyan plots the mean response. 
DOI: 10.7554/eLife.03427.013 

The following figure supplement is available for figure 6: 

Figure supplement 1. Numerical simulations of Dbxl neuron laser ablation experiments. 
DOI: 10.7554/eLife.03427.014 
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Sequentially deleting neurons in the model networks decreased frequency until the rhythm stopped 
altogether {Figure 6C,D, and Figure 6 — figure supplement 1B; Supplementary file 1). These 
simulations generally agreed with the experimental results except for two discrepancies. First, the 
amplitude of network output did not diminish (compare Figure 6C,E to Figures 2A and 3A). Second, 
rhythm cessation required the average deletion of 41 constituent model neurons (SD 1 5, SEM 6, range 
19-67, Figure 6 — figure supplement 1B,C', Supplementary file 1) as opposed to the experimental 
cell ablation tally of -85. First, we examine the loss of rhythmic function, and then address these 
discrepancies. 

To assess whether a collapse of network structure could explain the breakdown in rhythmicity, we 
computed canonical local and global measures of topology for the graph G(n,p) underlying each net- 
work simulation. These measures were computed after each cellular deletion and thus tracked contin- 
uously in parallel with the simulated networks. The table in Supplementary file 2 reports the value of 
each topological measure prior to any deletions and after the final deletion associated with rhythm 
cessation (see 'Materials and methods' for definitions and computational methods). From the start to 
finish, the cumulative ablation sequence caused no major change in local metrics including cluster 
coefficient, closeness centrality, and betweenness centrality. Global connectivity metrics such as the 
K-core, which has been applied to analyze rhythmic neural systems including the preBotC {Schwab 
et aL, 2010), showed only modest changes that were incommensurate with the large changes in 
frequency observed in experiments and simulations. The number of strongly connected components 
in the model networks did not depart from unity, thus the underlying graph was not fractured and 
every constituent interneuron could be reached via a finite number of synaptic links from every other 
interneuron, even after the rhythm stopped. These calculations show that these relatively low numbers 
of cumulative cellular ablations do not disconnect or disintegrate the core CPG, which suggests that a 
breakdown in network structure cannot explain the impairment and cessation of rhythmic function. 
The alternative is that neurons and synapses confer non-linear functional properties to the underlying 
rhythmic system that are not captured by the graph connectivity alone (see 'Discussion'). 

Dbx1 preBotC neurons with premotor function 

Lower cell ablation tallies perturbed and stopped the rhythm in simulations, and the aggregate burst 
magnitude did not decline {Figure 6C,E, and Figure 6 — figure supplement 1C, Supplementary file 1). 
Both disparities could be explained if a subset of the experimentally lesioned population consists of 
premotor — rather than rhythmogenic — interneurons. Thus, we tested whether Dbxl preBotC neurons 
project to the XII motor nucleus. Of the eight Dbxl neurons with inspiratory modulation {Figure 7A-D), 
two could be antidromically activated by stimulation within the XII nucleus. Figure 7 shows represen- 
tative data from such a Dbxl neuron whose Xll-evoked antidromic spike was extinguished by collision 
with an orthodromic spike triggered by a somatic current pulse {Figure 7E). Most Dbxl preBotC neu- 
rons are inspiratory and show commissural axons that cross the midline and innervate the contralateral 
preBotC {Figure BA-Q, as shown previously {Bouvier et aL, 2010] Picardo et aL, 2013). Here, we 
identify Dbx1 preBotC neurons that are also inspiratory modulated but send axons ipsilaterally toward 
the XII nucleus {Figure 8D-F and Figure 8 — figure supplement 1), consistent with a role related to 
premotor transmission of inspiratory drive from preBotC to XII motoneurons. 

Discussion 

Central pattern generators give rise to motor behaviors that are measurable in living animals as well 
as in reduced preparations that facilitate cellular-level investigations. To discover the core rhythmo- 
genic and premotor interneurons that comprise CPGs and establish their functional parameters remain 
formidable problems, particularly in mammalian systems. 

This report addresses the Dbxl core hypothesis, which posits that Dbxl preBotC interneurons con- 
stitute the rhythm-generating core for mammalian inspiratory rhythm. We affirm this hypothesis using 
cell-specific, cumulative laser ablation experiments, which also help quantify the size of the Dbxl 
preBotC population needed to defend inspiratory rhythmogenic function. Additionally, we provide 
unanticipated new evidence that Dbxl preBotC neurons also serve in a premotor capacity. While the 
existence of Dbxl premotor neurons in the preBotC is not a straightforward prediction of the Dbxl 
core hypothesis, it does not contradict it. We propose that an inspiratory rhythmogenic core popula- 
tion co-localizes with a subpopulation of premotor neurons, and both have the same developmental- 
genetic lineage related to Dbxl. These results help elucidate the rhythmogenic and premotor 



Wang et al. eLife 2014;3:e03427. DO!: 10.7554/eLife.03427 



11 of 25 



eLIFE 



Research article 



Neuroscience 



XII 



lOPioop preBotC 




lOPioop preBotC 




E 

-62 mV 
















XII stim 






app 



1 400 pA 



Figure 7. Dbxl preBotC neurons with premotor function. (A) Fluorescence and (B) bright field images of a slice 
preparation. Anatomical landmarks are illustrated including: XII, the hypoglossal motor nucleus; scNA, semi- 
compact division of the nucleus ambiguus; lOPioop, the ventral loop of the principal inferior olive, and the ventral 
surface of the preBotC. Scale bar is 1 00 pm and applies to A and B. A patch-recording pipette is visible, marking 
the inspiratory-modulated neuron detailed in C-E. A dotted circle indicates the tip of the pipette and cell body. 
(C and D) tdTomato expression, intracellular dialysis of Alexa 488, and merged image (C) from the inspiratory 
neuron shown with XII nerve output (D). Voltage and time calibration bars represent 20 mV and 1 s. Baseline 
membrane potential in the recorded neuron was -60 mV. (E) Antidromic activation of the Dbx1 inspiratory neuron 
from C and D. Action potentials were evoked by XII stimulation (left) and intracellular 5-ms supra-threshold current 
pulses (middle). When the antidromic XII stimulus was preceded immediately by a supra-threshold intracellular 
current pulse, the antidromic spike was occluded (collision test, right). Several sweeps, all from a -62 mV baseline 
membrane potential, are superimposed with vertical offset in each case. Voltage calibration is the same as panel D. 
Applied current (/app) calibration is shown. Time calibration bar for E is 25 ms. 
DOI: 10.7554/eLife.03427.015 
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Figure 8. Commissural and premotor projections of inspiratory Dbxl preBotC neurons. (A) Biocytin-filled and 
reconstructed Dbx1 preBotC neuron with commissural axon projection. The axon, which meanders in depth in this 
confocal image stack, was digitally traced (yellow) and superimposed in one plane for display. Axon trajectory 
crosses the midline of the slice and enters the preBotC contralaterally. Scale bar is 25 pm. (B) Mosaic image of the 
entire slice. The biocytin-filled soma (green) of neuron in A is shown at lower right (white arrow). Scale bar is 200 
[jm. Panels A and B have exactly the same orientation (dorsal up, ventral down). (C) Inspiratory discharge from the 
neuron in A and B. Top trace is membrane potential of the recorded Dbx1 preBotC neuron. Lower trace is XII 
output. Scale bars are 10 mV and 0.5 s. (D) Biocytin-filled and reconstructed Dbx1 preBotC neuron that projects 
toward the XII motor nucleus. Scale bar is 25 pm. The axon remained largely coplanar and thus is readily visible, 
except near its distal tip. In Figure 8 — figure supplement 1, this same neuron is shown with a digitally traced 
(yellow) axon superimposed on the confocal image. (E) Mosaic image of the entire slice. Neuron in D is shown at 
lower left (white arrow). Scale bar is 200 pm. Panels D and E have exactly the same orientation (dorsal up, ventral 
down). (F) Inspiratory discharge from the neuron in D and E. Top trace is membrane potential of the recorded Dbx1 
preBotC neuron. Lower trace is XII output. Scale bars are 10 mV and 0.5 s. 
DOI: 10.7554/eLife.03427.016 

The following figure supplement is available for figure 8: 

Figure supplement 1. Magnified view of the Dbxl preBotC neuron from Figure 8D-F in which the axon has been 
digitally traced in the confocal stack and superimposed over the image to better illustrate the axon projection 
toward the XII motor nucleus. 
DOI: 10.7554/eLife.03427.017 
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components of a key mannnnalian CPG, and furthermore provide quantitative parameters that govern 
its functionality. 

Photonic interrogation of CPG structure and function 

Two-photon lasers can destroy cells of a well-defined class with minimal damage to surrounding tissues 
{Eklof-Ljunggren et aL, 2012; Wang et aL, 2013). Here, we use this technique to study the contribu- 
tion of Dbx1 neurons in slices that capture essential components of the breathing CPG and generate 
measurable motor nerve output. Target detection relies on native fluorescent protein expression. An 
overwhelming majority of Dbx1 neurons in the ventral medulla have a glutamatergic transmitter 
phenotype and inspiratory modulated firing patterns {Bouvier et al., 2010] Gray et al., 2010] Picardo 
et aL, 2013), so the Cre/lox Dbx1 reporter system is a reliable means to identify neurons with inspira- 
tory function and target them for laser ablation. 

Dbx1 is also expressed in rostral parts of the ventral respiratory column between the caudal pole 
of the facial nucleus and the preBotC along the anterior-posterior axis {Feldman et aL, 2013] Gray, 
2013). The ventral respiratory column contains auxiliary inspiratory neurons {Figure 4) {Smith et aL, 
1990] Ballanyi et aL, 1999] Barnes et aL, 2007), which served as a control population. Laser ablating 
these neurons that do not have significant rhythmogenic function facilitates a comparative analysis of 
Dbx1 neuron ablations at the level of the preBotC. Laser ablation of Dbx1 neurons in the ventral res- 
piratory column had a minor effect on XII motor output amplitude and negligible effects on frequency 
and regularity. These negative results show that laser-tissue interactions are not generally deleterious 
for respiratory function in vitro {Eklof-Ljunggren et aL, 2012] Wang et aL, 2013). 

Dbx1 preBotC neurons form the inspiratory rhythmogenic core 

Periodicity is the hallmark feature of an oscillator. Here, sequential laser ablation of Dbx1 preBotC 
neurons steadily diminished the inspiratory burst frequency, caused cycle period fluctuations, and ulti- 
mately the cessation of rhythmic motor output. We conclude that the oscillator was continuously de- 
graded until it could no longer sustain spontaneous function. These data strengthen the proposal that 
Dbx1 neurons comprise the core inspiratory rhythm generator, which was originally based on Dbx1 
knockout mice that fail to breathe at birth {Pierani et aL, 2001), and an array of neuroanatomical and 
physiological criteria including glutamatergic transmitter phenotype, the expression of peptides and 
peptide receptors, strong inspiratory rhythmic phenotype, and the ability to synchronize the preBotC 
bilaterally {Bouvier et aL, 2010] Gray et aL, 2010] Picardo et aL, 2013). 

We previously laser-ablated rhythmic preBotC neurons identified by Ca^^ imaging. In that study, 
deleting all the detected targets (120 on average) slowed, destabilized, and then stopped the rhythm 
{Hayes et aL, 2012). The interpretability of these experiments suffered two caveats: the rhythm 
stopped after a delay of -30 min following the final target ablation, and furthermore, the transmitter 
phenotype of the ablated targets was unknown. In this study, use of the Dbx1 Cre-driver line 
ensured that the target neurons were glutamatergic, a requisite characteristic for respiratory 
rhythmogenic function (Greer et aL, 1991] Funk et aL, 1993] Ge and Feldman, 1998] Shao et aL, 
2003] Wallen-Mackenzie et aL, 2006). And here, destroying an average of 85 Dbx1 neurons stopped 
the XII motor rhythm in the midst of the ablation phase, before exhausting the target list, which 
suggests a more direct impact on the core rhythmogenic circuit. We cannot rule out the possibility 
that subsets of preBotC neurons generate 'burstlets' observable in local field recordings {Kam et aL, 
2013a). However, there is no collective inspiratory motor output after sequential laser ablation of 
Dbx1 preBotC interneurons, which indicates that the CPG is nonfunctional. Here, the debilitating 
effects on respiratory rhythm generation at a much lower ablation tally suggest that the preBotC core 
in vitro is very sensitive to the loss of just a few constituent interneurons (i.e., <100). This sensitivity to 
neuron loss may be accentuated in reduced slice preparations lacking excitatory and neuromodulatory 
drive from the rostral medulla and pons, as well as peripheral chemosensory and mechanosensory 
feedback via the vagus nerve. Extrinsic sources of drive raise preBotC network excitability and enhance 
respiratory rhythm. An acute sensitivity to neuron loss, such as we report for slices, may not apply to 
the preBotC network in vivo, but this remains to be tested via quantitative cellular ablation experiments 
with physiological monitoring in intact animal models. 

Parameters of the preBotC rhythmogenic core 

We detected an average of 705 Dbx1 target cells in preBotC-surface slices, but we conclude that a 
significant number were non-rhythmogenic. Some fraction of the detected targets can be discounted 
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as Dfox7-derived non-rhythmogenic glia {Gray et al., 2010). However, more significantly, some fraction 
manifests premotor function. The present evidence for premotor function in Dbx1 preBotC neurons with 
verified inspiratory discharge patterns (e.g., Figures 7 and 8D-F) is consistent with large-scale pressure 
ejections of biocytin in the preBotC region of another strain of Dbx7-reporter mice (Dbx7'-^^^ knock-in), 
which labeled many midline-crossing axons as well as axons projecting to the XII nucleus {Bouvier et aL, 
2010). Even 20 years ago it was recognized that a fraction of the excitatory neurons in the preBotC, and 
immediately dorsal to preBotC, had premotor functionality (Funk et aL, 1993). Because laser ablations 
in preBotC-surface slices decreased XII magnitude (e.g.. Figures 2A and 3A), we propose that a non- 
negligible fraction of the ablated Dbx1 neurons were inspiratory modulated but non-rhythmogenic, and 
most likely constitute XII premotor neurons [Peever et al., 2002) Chamberlin et aL, 2007] Koizumi 
et aL, 2008] Volgin et aL, 2008). This scenario explains why there was a decline in XII motor output in 
experiments (deletion of Dbx1 premotor neurons causes motor output to decline) that was not mimicked 
by model simulations of a pure rhythmogenic circuit. It also explains why sequential ablations in 
simulations perturbed and stopped the rhythm at much lower cell ablation tallies; a significant fraction of 
the photonically ablated Dbx1 neurons were unrelated to rhythmogenesis per se. 

Although we lack quantitative certainty, if we assume that each of the two caveats above (i.e., the 
existence of Dbx1 glia and premotor neurons) explains -10% of the detected targets, then the size of 
the essential preBotC core would be N = 705 - 2 [0.1 (705)] = 564, which is remarkably close to the 
estimate of -600 from adult rat studies that enumerated the population size based on NK1R expres- 
sion in the preBotC {Gray et aL, 2001, 1999). In our previous laser ablation study, we estimated pop- 
ulation size to be -325 {Hayes et aL, 2012), which probably underestimates the population size 
because incomplete fluorescent Ca^^ dye loading in slices precludes the detection of a significant 
fraction of the rhythmogenic preBotC network. 

Physiological significance 

Our data suggest that the preBotC contains rhythmogenic and premotor interneurons that both 
develop from Dbx7-expressing precursors. It is surprising that Dbxl neurons play these two roles in 
respiration when the role of Dbxl neurons in spinal locomotor systems seems to be coordinating left- 
right limb alternation at any speed {Lanuza et aL, 2004] Talpalar et aL, 2013) rather than rhythm 
generation or premotor transmission. In that regard, Dbxl preBotC neurons appear to have more in 
common with excitatory Shox2 interneurons of the lumbar spinal cord (a subset of V2 interneurons), 
which contribute to locomotor rhythm generation and premotor circuits downstream of the rhythm- 
generating core {Dougherty et aL, 2013). 

The present measurements imply that destroying on average X = 85 of N = 564 Dbxl preBotC 
neurons (15%) precludes spontaneous respiratory motor rhythm in vitro. The mean and its 95% 

confidence intervals are expressed as follows: X± Z^/2~T^ ' where is the cutoff value for a 

I Vn ) 

two-tailed normal distribution with probability a = 0.05, and is standard error. Thus, we conclude 

Vn 

the ability of Dbxl preBotC neurons to spontaneously generate rhythmic respiratory motor output in 
slice preparations is sensitive to deletion of 85 ± 1.96(20) Dbxl interneurons, i.e., 8-22% of its core. 
Although it cannot function spontaneously post-lesion, the rhythmogenic preBotC core does not 
appear to be destroyed by piecewise lesioning. Peptide injections evoked irregular transient bursts in 
lesioned preBotC-surface slices. Also, measures of local and global connectivity in lesioned network 
models remained undiminished after cumulative cell deletions stopped and precluded rhythmic func- 
tion. Therefore, we propose that preBotC function depends on non-trivial properties that emerge from 
non-linear synaptic and intrinsic membrane properties. Although such properties remain to be defini- 
tively determined, we advocate — and have explicitly modeled — a 'group pacemaker' rhythmogenic 
mechanism. In a group pacemaker, each constituent neuron forms recurrent connections with other 
constituent neurons in a network of finite size and amplifies excitatory drive via synaptically triggered 
inward currents {Rekling et aL, 1996] Rekling and Feldman, 1998] Rubin et aL, 2009). If that is a 
viable explanation for rhythmogenesis, then it could account for the loss of spontaneous function in 
the laser ablation context. Far before the cell ablation tally destroys the underlying network and its 
connectivity, the removal of each constituent neuron that contributes to rhythmic burst generation 
through its ability to amplify synaptic drive has a profound and deleterious effect on network functionality. 
The likelihood that arrhythmic respiratory networks retain considerable numbers of constituent neurons 
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and interconnectivity suggests that unraveling the cellular and synaptic mechanisms of rhythmogenesis 
and motor output could be exploited to restore functionality in lesioned slices and, to the extent that our 
observations apply in vivo, to develop clinical therapies that bolster respiratory function in pathological 
conditions of animal models or human patients. 

Dbx1 respiratory neurons in the medulla represent excellent potential targets for pharmacological 
intervention or gene therapy to treat respiratory pathologies. Potentially enhancing premotor func- 
tionality in Dbx7-derived neurons could ameliorate obstructive sleep apnea. Boosting the function of 
rhythmogenic Dbx1 neurons may mitigate central apneas of prematurity as well as opiate respiratory 
depression. Treatment strategies aimed at rhythmogenic Dbx1 neurons may help overcome the effects 
of a reduced quantity or efficacy of neurons within the preBotC due to neurodegenerative diseases or 
aging {Benarroch, 2003, Benarroch et aL, 2003, Tsuboi et aL, 2008). 

Materials and methods 
Ethical approval 

The Institutional Animal Care and Use Committee at The College of William & Mary, which ensures 
compliance with United States federal regulations concerning care and use of vertebrate animals in 
research, approved the following protocols. The anesthesia and surgery protocols are consistent with 
the 201 1 guidelines of the Animal Research Advisory Committee, which is part of the Office of Animal 
Care and Use of the National Institutes of Health (Bethesda, MD). 

Animal models 

We used transgenic mice that express Cre recombinase fused to the tamoxifen-sensitive estrogen 
receptor (CreERT2) in cells that express the Dbx1 gene {DbxV^^'^^^'^^) {Hirata et aL, 2009, Gray et aL, 
2010: Picardo et aL, 2013). DbxV^^''^^^'^^ mice were coupled to floxed reporter mice whose Rosa26 
locus was modified by targeted insertion of a /oxP-flanked STOP cassette followed by tandem dimer 
(td) Tomato {Gt(ROSA)26Sor^'°''-''°P-''^'^°'"''°, i.e., Rosa 26''^'^°"'''°, Jax No. 007905) {Madisen et aL, 2010). 
Tamoxifen administration to pregnant females on the tenth day after the plug date produces bright 
native fluorescence in Dbx7-derived neurons (i.e., Dbx1 neurons) in -50% of the offspring: DbxV^^'^^^^^', 
Rosa26^'^^°'^^^° . Dbx1 neurons can be visualized via native fluorescence in the preBotC and contiguous 
regions of the medulla. The DbxV^^'^^^^^ heterozygous line has a CD-1 background. The Rosa26^^^°'^^^° 
line is homozygous with C57BL/6J background. We verified animal genotype via real-time PCR using 
primers specific for Cre and tandem dimer red fluorescent protein. 

Neonatal mouse slice preparations 

Neonatal pups aged postnatal days 0-5 (PO-5) were anesthetized for at least 4 min of immersion in 
crushed ice in order to render the animals insentient to the same degree as would occur with inhala- 
tion anesthetics (Danneman and Mandrell, 1997; Fox et aL, 2007). Anesthesia via hypothermia 
facilitates the rapid isolation of the intact brainstem and spinal cord, which would otherwise be dam- 
aged by cervical dislocation. The brainstem and spinal cord were removed within 90 s and then dis- 
sected in a dish filled with artificial cerebrospinal fluid containing (in mM): 124 NaCI, 3 KCI, 1.5 CaCl2, 
1 MgS04, 25 NaHCOs, 0.5 NaH2P04, and 30 D-glucose, equilibrated with 95% O2 and 5% CO2 (pH = 7.4). 
After removing the meninges and isolating the XII nerve rootlets, the brainstem and contiguous 
upper cervical spinal cord were fixed in position on a paraffin-coated paddle, or glued to an agar 
block, with rostral side up. The paddle or block was mounted to the vise of a vibrating microtome. The 
advancing blade approached the ventral surface of the tissue for sectioning in the transverse plane. XII 
nerve rootlets remained visible during the sectioning sequence. We cut a single slice of thickness 
400-450 pm, which invariably retained the preBotC, XII premotor neurons and XII motoneurons that 
modulate and control airway resistance during breathing. 

We employed two discrete slice-cutting strategies to differentially expose respiratory nuclei at the 
slice surface, as previously described {Hayes et aL, 2012). The first slice type exposed the preBotC at 
the rostral face, and thus is called a preBotC-surface slice. The second slice type exposes the ventral 
respiratory column -100 jjm rostral to the preBotC at the rostral slice surface and served as a control 
slice for laser ablations. Histology atlases for newborn mice were used to calibrate slices online during 
sectioning {Ruangkittisakul et aL, 2011, 2014). For premotor recording experiments {Figure 7), we 
modified the preBotC-surface slice for the whole cell recordings such that the preBotC was exposed 
on the caudal surface. 
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Slices were perfused with 27°C ACSF at 4 ml/min in a recording chamber on a fixed stage upright 
nnicroscope. The external concentration was raised to 9 mM and inspiratory motor output was 
recorded from XII nerve roots using a suction electrode and an AC-coupled differential amplifier. The 
amplified electrical signal and a root-mean-squared (smoothed) version of the signal were recorded by 
a 16-bit analog-to-digital converter and stored on a digital computer. 

Because the composition of neural circuits at the rostral surface of the slice is critical for data 
interpretation, we fixed and stained each slice used for ablations at the end of the experiment to 
more precisely benchmark the neuroanatomical boundaries of respiratory-related nuclei according 
to the respiratory brainstem mouse atlases referred to above {Ruangkittisakul et a\., 2011, 2014). 
Fixation solution contained 4% paraformaldehyde in phosphate buffer (33 mM NaH2P04 and 67 
mM Na2HP04, pH = 7.2). After 1-hr in fixation solution, slices were rinsed in phosphate buffer for 
2 min, and then submerged for 60-75 s in staining solution containing 1% thionin acetate, 0.1 M 
sodium acetate trihydrate, and 0.1 M acetic acid. After washing in a series of ethanol solutions, 
slices were mounted in a well slide, obliquely illuminated, and digitally imaged via stereomicro- 
scope. Control slices were characterized by the compact division of the nucleus ambiguus (cNA), 
a thick dorsal inferior olive (lOD), and a minimally developed principal loop of the inferior olive 
(lOPioop). The preBotC-surface slices were characterized by very little (if any) visible portion of the 
cNA yet a clear semi-compact nucleus ambiguus (scNA), a fully developed lOPioop, as well as medial 
inferior olive (lOM). 

Electrophysiology 

We performed whole cell recordings using a Dagan (Minneapolis, MN) 1X2-700 current-clamp ampli- 
fier. Patch pipettes were fabricated from borosilicate glass (OD: 1.5 mm, ID: 0.87 mm, 4-6 MQ in 
the bath) and filled with solution containing (in mM): 140 K-gluconate, 10 HEPES, 5 NaCI, 1 MgCl2, 
0.1 EGTA, 2 Mg-ATP, 0.3 Na-GTP, 50-|jM Alexa 488 hydrazide, and 2-mg/ml biocytin. Empirical meas- 
urement of the liquid junction potential was 1 mV and thus not corrected. Access (series) resistance 
was -10-15 MQ, which was countered by bridge balance. Conventional current-clamp analog 
recordings were digitized at 4 kHz with a 1 6-bit A/D converter after 1 kHz low-pass filtering (PowerLab, 
AD Instruments, Colorado Springs, CO). 

Neurons were selected for recording based on native tdTomato fluorescence in neurons preferen- 
tially in the dorsal preBotC. After identifying an inspiratory Dbxl preBotC neuron, we tested for anti- 
dromic activation using a concentric bipolar electrode (FHC Inc., Bowdoin, ME) placed at the surface 
of the XII nucleus. Stimuli were triggered by a pulse generator (Tenma TGP1 1010 MHz Pulse Generator, 
Aim-TTi USA, Fairport, NY) and amplitude and polarity were controlled by a stimulus isolation unit 
(Iso-Flex, AMPI, Jerusalem, Israel). We applied cathodic stimuli at increasing intensities, to a maximum 
of 0.4 mA, to elicit short latency antidromic action potentials. Then, brief (1 ms) current pulses, at 
magnitudes at or exceeding rheobase, were applied before the antidromic stimulation such that both 
ortho- and antidromic spikes were evoked. The delay between stimuli was progressively decreased 
until a collision was observed, i.e., the antidromic spike was occluded. 

Biocytin cell reconstruction 

Biocytin-loaded neurons were fixed in 4% paraformaldehyde in 0.1 M Na-phosphate buffer for at 
least 16 hr at 4°C. Then, the slices were treated with Scale solution containing 4 M urea, 10% 
(mass/volume) glycerol and 0.1% (mass/volume) Triton X-100, for 10 days to clear the tissue and 
remove opaque background staining {Hama et a\., 2011). Slices were then washed in phosphate 
buffered saline (PBS) for 1 hr, followed by a 1 5-min cycle with PBS containing 1 0% heat-inactivated 
fetal bovine sera (F4135; Sigma-Aldrich). Next, slices were incubated in PBS containing fetal 
bovine sera with additional 1% Triton X-100. Finally, the slices were incubated in FITC (i.e., 
fluorescein-isothiocyanate)-conjugated ExtrAvidin (E2761; Sigma-Aldrich) overnight at 4°C, and 
then rinsed twice with PBS, followed by six 20-min washes in PBS, and then cover-slipped in 
Vectashield (H-1400 Hard Set, Vector Laboratories, Burlingame, CA). We visualized recorded 
neurons using a laser-scanning confocal microscope (Zeiss LSM 510, Thornwood, NY) or a spinning- 
disk confocal microscope (Olympus BX51 , Center Valley, PA). Images were contrast enhanced and 
pseudo-colored using the free /mageJ software (National Institutes of Health, Bethesda, MD), and 
then digitally reconstructed using the free Neuromant/c software for morphological reconstruction 
(Myatt et a\., 2012). 
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Laser ablation: target detection 

Dbx1 neurons were detected and mapped within three-dimensional (3D) volumes of the preBotC or 
ventral respiratory column, and then subsequently laser ablated while monitoring respiratory network 
functionality. The instrument incorporated a Zeiss LSM 510 laser scanning head and fixed-stage 
microscope body with a 20x/1.0 numerical aperture water-immersion objective, an adjustable 
wavelength 1.5 W Ti:sapphire tunable laser (Spectra Physics, Irvine, CA), and a robotic xy translation 
stage (Siskiyou Design, Grants Pass OR). The methodology has been described in a technical report 
(Wang et aL, 2013) and in an original research report {Hayes et aL, 2012). 

We wrote custom software dubbed Ab/ator that automated a three-step routine. The first step 
(initialization phase) defines the domain for target detection and ablation. The domain can be bilaterally 
distributed, like the preBotC and ventral respiratory column. The maximum size of any part of the domain 
in the transverse (xy) plane must fit within an area of maximum dimensions 412 square micrometers. 
The z domain (depth) is a function of tissue opacity, laser power (Tiisapphire), and the emission properties 
of the fluorescent reporter. For neonatal mouse brainstem tissue (PO-5), using 800-nm pulses emitted at 
~1 W, which measured 36 mW at the specimen plane, the z domain generally measured less than 1 00 pm. 

The second step (detection phase) acquires high-resolution images via confocal microscopy with a 
visible-wavelength laser (HeNe 543 nm for tdTomato). Dbxl neurons were identified by native fluorescent 
protein expression using a threshold-crossing target detection algorithm in Ablator software, which is 
open-source and available for free download at the sourceforge.net archive, i.e., http://sourceforge.net/ 
projects/ablator/. Additional image processing routines differentiate Dbxl somata from auto-fluorescent 
debris and neuropil (Figure 1 — figure supplement 2). The final map of Dbxl neuron targets reflects the 
position of the center of each cell body in the 3D volume of the domain (see Figure ID). 

Iterative threshold-crossing algorithm and innage processing 

Given an image that captures features of potential targets. Ablator calls the Analyze Particles routine 
in ImageJ, which is free image analysis software in the public domain {Schneider et al., 2012), to 
select ROIs. This routine detects particles using a threshold based on pixel intensity. It starts with a 
high value (near the maximum) and iteratively drops the threshold while accumulating particles (i.e., 
target ROIs) within a certain range of areas specified within the Ablator configuration. At first, with a 
high threshold, few local maxima are detected and the mask is small and sparse with ROIs. The routine 
then lowers the threshold by a user-defined increment and re-analyzes the image. As threshold decreases 
in steps, more ROIs become detectable. These newly detected ROIs are added to the mask, which 
expands the list of potential targets. The threshold is decreased incrementally over a number of partitions 
determined by the quotient of 4096 values of fluorescence intensity (for a 12-bit image) divided by 
the user-defined increment (above). As the detection process continues and threshold decreases incre- 
mentally, the smaller ROIs from prior iterations — which are fully contained in larger ROIs from the current 
iteration — are discarded, and the new larger ROIs are retained. Conversely, if a newly detected ROI at 
the current iteration envelopes two or more ROIs from a prior iteration, then the newly detected superset 
ROI is discarded and the multiple ROIs from the earlier iteration are retained. Thus the system avoids 
spuriously conflating two (or more) cells into a single target. After looping through all the partitions, the 
remaining set of ROIs is saved as the mask of potential target neurons for that focal plane. 

Circularity test 

Ablator evaluates the circularity of ROIs as part of the threshold-crossing algorithm. A circularity score 

C = 4/7"-^, is computed based on a (area) and p (perimeter) of the ROI. a and p are measured by the 
P 

Analyze Particles routine in ImageJ {Schneider et al,, 2012). C ranges from 0 to 1. Scores near 0 
denote an elongated polygon. C approaches 1 for a perfect circle. If C falls below a user-specified cut- 
off, then the ROI is rejected from the target list. The appropriate C score depends on the characteristic 
morphology of the neurons of interest. Valid Dbxl neurons in the preBotC and ventral respiratory 
column pass the circularity test when C exceeds 0.75. Circularity is particularly useful in selecting 
somata rather than neuropil or auto-fluorescent detritus as valid targets {Figure 1 — figure supplement 2). 
Rejecting an isolated dendrite segment {Figure 1 — figure supplement 2B,Q avoids protracted lesion 
attempts during the ablation phase of the experiment, which are problematic because the dendrite 
and its soma are redundant targets and attempting to laser-ablate the dendrite is more likely to sever 
the process rather than kill the neuron {Kole, 2011). 
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Priority rule 

This final processing step eliminates redundant targets from adjacent focal planes. The same cell tar- 
gets can be detected, and pass the circularity test, in more than one plane. When overlaying ROIs exist 
within adjacent focal planes the earliest acquisition from the deeper plane is retained and all other 
superficial ROIs are deleted {Figure 1 — figure supplement 2C,D). 

Cell-specific laser ablation 

Ablator chooses Dbx1 neuron targets in random order and advances until all the targets are exhausted 
or the respiratory rhythm ceases for longer than 120 s. The Tiisapphire laser scans a 10 square micro- 
meter spot centered on each target with 800-nm pulses at maximum intensity. The ablation is con- 
firmed if fluorescence is detected in the band 560-615 nm, which reflects presumed water vapor in the 
cell cavity and excludes infrared reflections of the long-wavelength laser {Figure 2 — figure supple- 
ment 1A) {Wang et al., 2013). In addition, lesioned targets disappear from the fluorescence image 
{Figure 2 — figure supplement IB), and their pre-lesion bright field image {Figure 2 — figure supple- 
ment 1C) is replaced post-lesion by a pock mark {Figure 2 — figure supplement ID). Confirmed 
lesions add to a running tally. If lesion confirmation cannot be obtained, then the target selection 
algorithm does not advance and subsequent attempts are made to lesion the ROI. With each subse- 
quent iteration, the scanning speed is decreased to improve the likelihood of lesioning the target. This 
loop repeats a total of five times. If confirmation of lesion cannot be ascertained after the fifth attempt, 
then it is deemed a failed lesion. Failed lesions do not contribute to the tally and their ROIs are 
removed from the list of targets to avoid reselection for the remainder of the experiment. A log file 
documents lesions by index number and time of confirmation. The XII rhythm is monitored and 
recorded continuously so its state can be directly correlated with the lesion tally in real time. Cell tar- 
gets are destroyed in successful lesions so their effects are cumulative. The laser lesions are performed 
bilaterally in the preBotC. After a batch of lesions on one side, the robotic xy translation stage translates 
to the contralateral side and performs another batch, and then switches sides again, and so on until 
the targets are exhausted or the XII rhythm ceases. 

We measured XII burst magnitude (amplitude and area) and computed cycle period (the interval 
between consecutive XII bursts) using LabChart software (ADInstruments, Colorado Springs, CO). The 
regularity score (RS) was defined as the quotient of period of the present cycle Tn with respect to the 
mean cycle period for ten previous cycles: 

ly^ T • 

j Z./=1 ' n-i 

where j = 10. We defined the control epoch as 30 min of continuous recording from the end of 
the detection phase to the beginning of the ablation phase. Data sets were tested for normality using 
a Shapiro-Wilk test. We rejected the null hypothesis that the data are drawn from a normal distribution 
if the p-value of the test statistic was less than a = 0.05. Data that could be considered normally distrib- 
uted were compared using two-tailed paired t-tests, whereas data that did not conform to the normal 
distribution were compared using non-directional (two-tailed) Mann-Whitney (J-tests. XII burst 
amplitude and frequency/cycle period were reported with standard deviation (SD) and standard error 
of the mean (SEM). Discrete cell counts that pertain to the number of neurons detected or the number 
of neurons lesioned are reported with SD, SEM and min-max range. 

Network simulations and modeling 

We wrote a Matlab (MathWorks Inc., Natick, MA) script to generate Erdos-Renyi G(n,p)-directed random 
graphs {Newman eta/., 2006) with key parameters of population size (n) and connection probability (p). 
Vertices (a.k.a., nodes) of G(n,p) were populated by Rubin-Hayes preBotC neuron models and the 
directed edges (a.k.a., links) between vertices were modeled by excitatory glutamatergic synapses 
{Rubin et al., 2009). We simulated the network models on the SciClone computing complex at The 
College of William & Mary, which features 193 nodes with a total of 943 CPU (central processing unit) 
cores, 5.9 terabytes of physical memory, 220 terabytes of disk capacity, and peak performance of 21.2 
teraflops. We used a Runge-Kutta fourth-order numerical integration routine with fixed time step of 0.25 
ms. Network models were subject to 100 random deletions, one deletion every 25 s. Neuron deletions 
were achieved by setting the synaptic state variable and its corresponding differential equation to zero. 
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which essentially removes the cell from the network. Deleted neurons no longer contributed to 
running-time histograms of network activity and were removed from raster plots (e.g., Figure 6C). 
Transient glutamatergic stimulation of constituent model neurons mimicked the experimental gluta- 
mate un-caging protocol by Kam et al. which evoked respiratory bursts in the preBotC {Kam et al., 
2013b). Focal stimulation was achieved by setting the synaptic state variable to 0.9 for 200 ms, 
without modifying the differential equation, so the glutamatergic excitation was indeed transient. 
Focal stimulation was applied to rhythmically active networks several seconds following an endoge- 
nous burst {Figure 6B). 

Since there is uncertainty regarding the exact network size, we conducted a series of simulations for 
a range of (n,p) with the aim of finding a reasonable parameter range to produce respiratory-like 
rhythms (3-4 s cycle period prior to ablations). We varied n from 200 to 400 with a step size of 10 and 
p from 0.1 to 0.2 with a step size of 0.0125. For each parameter set, 10 simulations without deletion 
were conducted for 25 s to assess network rhythmicity {Figure 6 — figure supplement 1). 

For the parameter sets whose initial period fell between 3 and 4 s, we performed 5-6 simulations 
with deletions (for n = 320, 330, 340 we performed 16 simulations in each case) and then calculated 
the longest period, the ablation tally, and discrete network metrics pertaining to G(n,p). The results are 
documented in Figure 6 and Figure 6 — figure supplement 1, as well as the table in Supplementary 
file 1. During the simulations, raster graphs were simultaneously generated to detect the spiking for 
each individual neuron {Figure 6Q. The running time histogram is based on the raster graph for each 
simulation, from which we computed the cycle period and amplitude (number of spikes per time bin. 
Figure 6D,E}. 

Discrete network simulations 

A network with n vertices can be represented by its adjacency matrix A(n x n) in a manner that if there 
is a connection from vertex / to vertex j then Ajj = 1, otherwise A,y = 0. The adjacency matrices are 
asymmetric for neuronal networks, which are directed (i.e., the chemical synapses are unidirectional). 
In discrete simulations, the lesion of neurons is modeled by removing vertices from the adjacency 
matrix along with their edges, i.e., connections (in and out). We computed three global metrics (K-core, 
number of strongly connected components, average in and out degree) for the initial network and the 
remaining network after a sequence of 100 random deletions. Also for each deleted vertex, we 
computed three local network metrics (local cluster coefficient, closeness centrality, betweenness cen- 
trality) to indicate the importance of the vertex within the previous network. The metrics are defined 
below and reported in the table of Supplementary file 2. 

K-core 

It refers to the maximum sub-graph such that each vertex of the sub-graph has at least K edges 
(connections). In this case, an in-arc and an out-arc both count as an edge. 

Strongly connected connponents (SCC) 

The strongly connected components of a directed graph G(n,p) are its maximal strongly connected 
sub-graphs, such that within each sub-graph there is a path from each vertex to every other vertex. 
Therefore the number of SCC can exceed unity. Nonetheless, when SCC = 1 the existing network is 
said to be fully connected, i.e., there are no isolated islands and every vertex can connect to every 
other vertex via a finite number of edges. 

Average in and out degree 

In an n X n adjacency matrix A of a directed graph G(n,p), Ajj = 1 refers to a connection from vertex / 

n 

to j, and A^, =1 refers to a connection from vertex jto /. Therefore y^A. is the out-degree for node 

n i=1 

/ while /A. is the in-degree. 

Local cluster coefficient 

Local cluster coefficient measures how close the neighbors of the vertex are to being a complete graph, 
i.e., a graph where each vertex is connected to every other vertex. For a vertex v, with k-, edges, the 
local cluster coefficient is defined as 

^ _ \{A.,■.v.,v,^N,A., = ^}\ 

/c,.(/c,.-1) 
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where N, is the neighborhood of v-,, the sub-graph formed by all the vertices v, connects to (that is, all 
the out-neighbors of v-,). The numerator is the number of actual connections within N-, while the 
denominator is the number of connections if N, is a complete graph. 

Closeness centrality 

For a vertex v-,, the farness is defined to be the sum of shortest paths from v-, to every other reachable 
vertex. The closeness of v-, is the inverse of the farness. The closeness centrality for v-, is defined as the 
product of the number of vertices in the graph n and the closeness of v,-. From this definition, a central 
vertex would have a small farness and a large closeness centrality. 

Betweenness centrality 

Betweenness centrality measures the frequency that a vertex acts as a bridge in the shortest path 
between two other vertices. It is defined as Cg (v) = ^ ^stM ^[^g^g v, t are three different vertices 

in the graph, and cr^tM is the number of shortest paths between s and t through v, while a^^ is the total 
number of shortest paths between s and t. Betweenness centrality is usually normalized by dividing 
the number of total possible vertex pairs (n - 1 )(n - 2), excluding v. 
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Additional files 

Supplementary files 

• Supplementary file 1. Numerical simulations of Dbx1 neurons in model preBotC networks subjected 
to cumulative laser ablation experiments. Erdos-Renyi random directed graphs G(n,p) were populated 
with Rubin-Hayes preBotC neuron models at each node, and their links were described by excitatory 
synapses, as described for Figure 6 and Figure 6 — figure supplement 1 above. In numerical simulations, 
the resulting network models with very high probability of generating rhythm and respiratory-like 
cycle period (~4 s, indicated by asterisks in Figure 6A and Figure 6 — figure supplement 1) were 
subjected to piecewise cumulative ablation protocols like slice experiments (Figures 2-4). The 
parameters describing the model networks (number of neurons n and synaptic connection probability 
p) are listed below in columns 1 and 2. Each ablation experiment was simulated five or more times. The 
maximum period (in s) and cumulative ablation tally (unitless) required to stop the rhythm are listed in 
the table for each individual realization of the network model along with average values for these 
characteristic measures. The networks deemed to be most representative of the preBotC, i.e., (n,p) = 
(320, 0.1375), (330, 0.125), (340, 0.125) were simulated 16 times each. 

DOI: 10.7554/eLife.03427.018 

• Supplementary file 2. Discrete network simulations. As in the main text and the table above, 
parameters (n,p) represent the number of constituent neurons and connection probability. Here 
the networks are Erdos-Renyi static directed random graphs G(n,p) (as in Figure 6 and Figure 6 — 
supplement 1); that is, the nodes are not populated with dynamical models and the interconnections 
between nodes are simply static directed links (rather than dynamical synapses). Each network (static 
graph) below was subjected to 100 random deletions. During the ablation sequence, we computed 
the following global network metrics: K-core, the number of strongly connected components (SCC), 
and local network metrics: local cluster coefficient, closeness centrality, and betweenness centrality. 
The initial (first) and final (last) values for the global and local measures are plotted side by side in the 
appropriate columns below. We also computed the initial in- and out-degree (i.e., the average number 
of directed connections in and out), percentage drop in the final average in-degree, and the 
percentage-drop of the final average out-degree. The 'Ave.' row reports average change (in percent) 
for K-core, SCC, cluster coefficient, closeness centrality, and betweenness centrality, as well as the 
average in- and out-degree for initial and final states of the network. Definitions for the characteristic 
measures are elaborated in the 'Materials and methods'. 

DOI: 10.7554/eLife.03427.019 
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